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Abstract: We present a three-loop model of neutrino mass in which both the weak scale 
and neutrino mass arise as radiative effects. In this approach, the scales for electroweak 
symmetry breaking, dark matter, and the exotics responsible for neutrino mass, are related 
due to an underlying scale-invariance. This motivates the otherwise-independent O(TeV) 
exotic masses usually found in three-loop models of neutrino mass. We demonstrate the 
existence of viable parameter space and show that the model can be probed at colliders, 
precision experiments, and dark matter direct-detection experiments. 
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1 Introduction 

Radiative symmetry breaking [1] offers an interesting alternative to the conventional Higgs 
mechanism. In this approach, calculable weakly-coupled radiative effects induce symmetry 
breaking in classically scale-invariant theories, thereby giving birth to mass — a process 
known as dimensional transmutation. When applied to the Standard Model (SM), it is well 
known that radiative symmetry breaking is not viable, due to the destabilizing influence of 
the heavy top quark. However, the SM is known to be incomplete, due to e.g. an absence of 
massive neutrinos and the need to incorporate dark matter (DM). It is therefore interesting 
to consider the viability of radiative symmetry breaking within SM extensions. 
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The addition of massive neutrinos and DM to the SM likely requires new degrees of 
freedom. When considering radiative symmetry breaking, there are a number of relevant 
considerations that can guide the choice of beyond-SM fields. The destabilizing radiative 
corrections from the top quark can be overcome by bosonic degrees of freedom with mass 
> 200 GeV. In principle these states could be much heavier than the TeV scale. However, ra¬ 
diative symmetry breaking typically introduces a single scale into a theory, with other mass 
and symmetry breaking scales related to this scale. ^ Consequently both the electroweak 
scale and the mass scale for exotics may be related via dimensionless parameters. Thus, 
absent hierarchically small parameters [2], one anticipates exotics with O(TeV) masses. 

In the LHC era, TeV scale exotics are of particular interest. However, efforts to generate 
tiny neutrino masses via weak-scale exotics can struggle to achieve the necessary mass- 
suppression, relative to the weak scale, without invoking tiny couplings. Perhaps the most 
obvious exception are models with radiative neutrino mass, as the inherent loop-suppression 
in such models can motivate lighter new physics. From this perspective, three-loop models 
of neutrino mass are particularly compelling, as the new physics is expected to be O(TeV). 

These considerations focus our attention on scale-invariant models with three-loop neu¬ 
trino mass. If we also seek to address the DM problem, a minimal approach would see the 
DM play a role in either generating neutrino mass or triggering electroweak symmetry 
breaking. Thus, we arrive at a picture in which both the weak seale and neutrino mass arise 
as radiative effeets, with the weak scale, the DM mass, and the mass scale for the exotics 
that induce neutrino mass, all finding a common birth, via dimensional transmutation. This 
picture can address short-comings of the SM, while also explaining why the exotics required 
in three-loop neutrino mass models have (otherwise independent) masses of O(TeV) — a 
common ancestry requires that they be related to the weak scale. 

In this work we present a scale-invariant model for three-loop neutrino mass that con¬ 
tains a fermionic DM candidate. We explore the model in detail and present feasible 
parameter space that achieves the correct DM relic abundance, while generating viable 
symmetry breaking and neutrino masses — all compatible with low-energy constraints. As 
per usual for scale-invariant frameworks, the model predicts a dilaton. However, here the 
dilaton has the dual role of allowing electroweak symmetry breaking and simultaneously 
sourcing the lepton number violation that allows radiative neutrino masses. We note that a 
number of earlier works studied relationships between the origin of neutrino mass and DM, 
see e.g. Refs. [3-8]. There has also been much interest in scale-invariant models in recent 
years, see e.g. Refs. [9-14]. 

The structure of this paper is as follows. In Section 2 we introduce the model and 
detail the symmetry breaking sector. We turn our attention to the origin of neutrino mass 
in Section 3 and discuss various constraints in Section 4. Dark matter is discussed in 
Section 5 and our main analysis and results appear in Section 6. Conclusions are drawn in 
Section 7. 


^The exceptions being when a theory also contains a confining gauge sector, as with QCD in the SM, or 
a completely decoupled hidden sector possessing its own symmetry breaking and/or confining pattern. 
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2 A Scale-Invariant Three-Loop Model 


We consider a classically scale-invariant (SI) extension of the SM in which neutrino mass 
appears at the three-loop level. The SM is extended by the addition of two charged 
scalars, 3^2 ~ (1)1)2), three singlet fermions, Aj/j ~ (1,1,0), with i S {1, 2, 3} la¬ 
beling generations, and a singlet scalar, (p ~ (1,1,0).^ A Z2 symmetry with action 
{52, Nr} —)• { — 52, —Nr} is imposed, with all other helds being Z 2 -even. This symmetry 
remains exact in the full theory, making the lightest Z 2 -odd held a stable DM candidate, 
which should be taken as the lightest fermion, Ni = A^M) to avoid a cosmologically-excluded 
stable charged particle. The scalar cp plays a key role in triggering electroweak symmetry 
breaking, as explained below, and also ensures that lepton number symmetry is explicitly 
broken, thereby allowing radiative neutrino mass. 

Consistent with the SI and Z2 symmetries, the Lagrangian contains the following terms: 

C^-{fo, 0L%LfiSt+giaMStec.R + ^-^} - \m(p'WfNi - A(A,5i,2,</-), (2.1) 

where Greek letters label SM havors, a, /3 S {e, //, r}, and fa/3, g/a and y/ are Yukawa 
couplings. The Z2 symmetry forbids the term LHNr, which would otherwise generate 
tree-level neutrino masses after the SM scalar H ~ (1, 2,1) develops a VEV. The potential 
V{H, S\^2,(p) is the most-general potential consistent with the SI and Z2 symmetries. 

2.1 Symmetry Breaking 

We are interested in parameter space where both (p and H acquire nonzero vacuum ex¬ 
pectation values (VEVs), (H) 0 and ((/>) 0. This breaks both the SI and electroweak 

symmetries while preserving the Z2 symmetry. The most-general scalar potential includes 
the terms 

Eo(A, 5 i, 2 , </>) D Ah + + + \{>^^a\H\^ + \^a(p')\Sa? 

a=l ,2 ^ 

( 2 . 2 ) 

A complete analysis of the potential requires the inclusion of the leading-order radiative 
corrections. In general the full one-loop corrected potential is not analytically tractable. 
However, a useful approach for approximating the ground state in SI models was presented 
in Ref. [15]. Taking guidance from Ref. [15], we adopt an approximation for the ground state 
that allows one to obtain simple analytic expressions. The physical spectrum contains two 
charged scalars 5]^2) two neutral scalars, denoted as hi^2- As discussed in Appendix A, 
for the present model, the minimum of the loop-corrected potential can be approximated by 
neglecting loop corrections involving only the scalars hi ^2- The viability of this simplihcation 
follows from the dominance of the beyond-SM scalars S}~2 (see Appendix A). Adopting this 

^Quantities in parentheses refer to quantum numbers under the SM gauge symmetry SU (3)c 0 SU ( 2 ) l 0 
U{1)y. 
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approximation, the one-loop corrected potential for the CP-even neutral scalars is 


Vi-l {h, (t>) 


G(r/) 


i=all fields 


GdTT^ 



(2.3) 

(2.4) 


where A is the renormalization scale, are the field multiplicities, and we employ the 
unitary gauge, with H = (0,/i/\/2)^- The sum is over all fields, neglecting the light SM 
fermions (all but the top quark) and the (to be determined) neutral scalar mass-eigenstates 
hi^ 2 - Due to the SI symmetry, the held-dependent masses can be written as 

'miiKcj)) = + (2.5) 

the constants Oj and /?* are given by 

+ g'"^ 2 

= Y’ oiz = -^-, at = yt, as^ = Xna, o^Ni = 0, 

Pw = Pz = Pt = 0, Psa = X(j,a, PNi = 2yf, (2.6) 


with g {g') and yt are the SU{2)l (I7(l)y) gauge and top Yukawa couplings, respectively. 

Dimensional transmutation introduces a dimensionful parameter into the theory in 
exchange for one of the dimensionless couplings. In the present model, an analysis of the 
potential shows that a minimum with {h) = u 7^ 0 and {p) = x 7^ 0 exists for A^h < 0, and 
is triggered at the scale where the couplings satisfy the relation 


2 


AhA,* H- TT 


PE 


rii { P; 




riiaiG' (mf) 
i 


0, (2.7) 


with G' (g) = dG (g) /dg. The further condition 


2 Ah 


+ Eir^G'K), 


Ah x'^ 


( 2 . 8 ) 


is also satisfied at the minimum. Thus, for A,^h,h = one has u ~ x and the exotic scale 
is naively expected around the TeV scale. Note that Eqs. (2.7) and (2.8) ensure that the 
tadpoles vanish.^ 

Defining the one-loop quartic couplings as 

g-z .1-1 .1-1 d^Vi-i 

^ 6 6 dh^ ’ dK^dcp^' ^ ’ 

vacuum stability at one-loop requires that the following conditions be satisfied: 

Ah"', ai-'. Ay + 2^\]r^\Y > 0. (2.10) 

^To our level of approximation, Eqs. (2.7) and (2.8) are the loop-corrected generalizations of the standard 
tree-level results, 4 ^Ah(A) \^{X) -|- A 0 h(A) = 0 and A<^h/ 2 Ah = v'^/x^ [16]. 
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We must also impose the condition < 0 to ensure that the vacuum with v ^ 0 and 
x 7^ 0 is the ground state. ^ Eq. (2.10) also guarantees that the eigenmasses-squared for the 
CP-even neutral scalars are strictly positive, and forces one of the beyond-SM scalars £'^*'2 
be the heaviest particle in the spectrum. 


2.2 The Scalar Spectrum 

The mass matrix for the neutral scalars is denoted as 

U_, 0 . 0) 3 (2,11) 

2 ) K'l’) 

where the mass parameters rrihh, rrifi)^ and ruhtf) are calculated from the loop-corrected 
potential Vi-i {h,(f)). The mass eigenstates are labeled as 


hi = cos Ofih — sin 9h4> ^ /12 = sin 9hh + cos 9h 4> j 

with the eigenvalues and mixing angles given by 


( 2 . 12 ) 


^4^,. = o 


tan20/i = 


mh + m^r, ± 


2771^2 


'-22 


77722 — m 


11 


f + ^rn\^ 


(2.13) 

77722 - "^11 

Here hi is a massive SM-like scalar and /i2 is a pseudo-Goldstone boson associated with SI 
symmetry breaking — the latter is massless at tree-level but acquires mass at the loop-level. 
One can obtain simple tree-level expressions for the SM-like scalar mass 

Ml = (2Ah - A^H)^^ (2.14) 


and the mixing angle, 

X V 

Ch = cos 9h = . , Sh = sm9h = —j===, (2.15) 

V a; + V V a: -|- 

though in large regions of parameter space it is important to include loop corrections to 
these expressions to obtain accurate results. In our numerical analysis we employ the 
full loop-corrected expressions for the scalar masses and mixing, as is necessary to obtain 
M/12 7^ 0- SI symmetry, the parameters in the model are somewhat constrained, 

with A^ and A^h fixed by Eqs. (2.7) and (2.8) while the Higgs mass M/^ ~ 125 GeV fixes 
Ah. 

The tree-level masses for the charged scalars, 3^2^ 

\ + Ana'i^^} for a = 1, 2, (2.16) 

where Sf and do not mix due to the Z 2 symmetry. Note that a useful approximation 
for M/12 f® [f^l 


< =* 8,7m‘+W0 ( K + 6Vh-- + + 2^M|, - 2^ , (2.17) 


a=l 


i=l 


'^For AJjjj > 0 the vacuum with only one nonzero VEV is preferred. 
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which shows that one of the beyond-SM scalars 51^2 must be the heaviest beyond-SM state 
in order to ensure > 0. 

As mentioned already, we expect the VEVs to be of a similar scale, {(j)) ~ (/i), as 
evidenced by Eq. (2.8). For completeness, however, we note that there is a technically 
natural limit in which one obtains (</>) S> {h). This arises when all the couplings to cf) are 
taken to be hierarchically small, namely {yj, A^hj ^ 01 , 2 } ^ Ij with the masses Mf^, 
and Ms ^2 held at C)(TeV). This feature reflects the fact that (f) decouples in the limit 
{yi) •^ 01 , 2 } 0, up to gravitational effects [17]. In this limit we expect the model to 

be very similar to the KNT model [3], but with a light, very weakly-coupled scalar in the 
spectrum, / 12 . Absent a compelling motivation for such hierarchically small parameters, we 
restrict our attention to values of {(j)) <5 TeV. 

3 Neutrino Mass 

We now turn to the origin of neutrino mass. The Z 2 -odd fermions, Ni, develop masses 
M^r. = yi{4>), and do not mix with SM leptons due to the Z 2 symmetry. We order their 
masses as Mqm = < ^In 2 < SM neutrinos, on the other hand, acquire mass 

radiatively. The combination of the Yukawa interactions in Eq. (2.1) and the term 

V{H, 5i,2, </>) D ^(5^)'(^2+)^ (3-1) 

in the scalar potential, explicitly break lepton number symmetry. Consequently neutrino 
masses appear at the three-loop level as shown in Figure 1. 

Si S 2 / \S2"^^Si 

' 

y' ! \ 

\ ^ 

-► f C ^ ^ ^4^ f ^^^—>■ ^ ^ 

el/ el Tr ^ Tl vl 

/ \ 

>7 X V 

W ( 0 ) W 


Figure 1. Three-loop diagram for neutrino mass in a scale-invariant model. 


Calculating the loop diagram, the mass matrix has the form 


(-Mj/)q^ 


As 

(4yr2)3 


m„mp 

Ms2 


* * i? 

9(ji dpi 11 


3p ^ P^loop 


w 


(3.2) 


where mo-,p denote charged lepton masses and the function Fioop{x,y) encodes the loop 
integrals [5] 


Floopip^i P') 


^ r dr ^ 

8/12 io ^ + 


, x{l — x)r + — x)P + X 

ax In - 


x{l — x)r + X 


(3.3) 
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One can relate the neutrino mass matrix to the elements of the Pontecorvo-Maki- 
Nakawaga-Sakata (PMNS) mixing matrix [18] elements, we parameterize the latter as 


U, 


f C12C13 C13S12 5136“*'^'^ ^ 

—C23'S12 — Cl2Sl3'S23e*‘^‘^ C12C23 — Si2Si3S23e*'^'^ C13S23 

\ S12S23 - Cl2C23Sl3e*^‘' -C12S23 - C23Sl2Sl3e*^‘' C13C23 / 


y. Urt 


(3.4) 


with 6d the Dirac phase and Um = diag(l, encoding the Majorana phase de¬ 

pendence. The shorthand Sij = sin0jj and Cij = cos Oij refers to the mixing angles. For our 
numerical scans (discussed below) we fit to the best-fit experimental values for the mixing 
angles and mass-squared differences: 5^3 = 0.3201qq}7, S 23 = 0.431qq 3, sfg = 0.0251 qqq 3, 
|Am 33 | = 2.55lQQg X 10“^eV^ and = 7.62lQ|g x 10“®eV^ [19]. Furthermore, we 

require that the contribution to neutrino-less double beta decay in this model satisfies the 
current bound. Within these ranges, one determines the parameter space where viable 
neutrino masses and mixing occur in the model. 


4 Experimental Constraints 

In this section we discuss the constraints on the model from the lepton flavor violating 
process /x —>• eq, the electroweak precision tests, the invisible Higgs decay, and the effect on 
/177 process. 


4.1 Lepton flavor 

Flavor changing processes like /x ^ e -|- 7 arise via loop diagrams containing virtual charged 
scalars and give important constraints on the model. At one-loop the branching ratio for 
/X —^ e -|- 7 is 


B{n e'j) = 


r(/x^e-t7) 
r(/x e + u + i>) 


av 

3847r 




IJ-T J Te\ 

Ml, 


+ 


36 


E 


glg,^F2{MflMl) 


(4.1) 


where ^ 2 ( 1 ?) = [1 — 6ii-|-3i?^-|-212^ — 612^ log i?]/[ 6 (l — i?)^]. The corresponding expression 
for Bt —)• /X -|- 7 ) follows from a simple change of flavor labels in Eq. (4.1). Similarly, the 
one-loop contributions to the anomalous magnetic moment of the muon are 


6a^ = 







6 M|^ 



Ml 


F^iMflMl 


(4.2) 


Null-results from searches for neutrino-less double-beta decay give an additional constraint 
of {Aiu)ee ^ 0.35 eV [20], though we find this is easily satisfied. 
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4.2 Electroweak precision tests 


In principle, precision electroweak measurements can provide additional constraints. The 
oblique parameters characterizing new physics effects are given by [ 21 ] 


a Azz (A'/f) - Azz (0) 

4s2 ^ ~ 


_ Aww ( 0 ) _ Azz ( 0 ) 
M| 


dAjj (q^) 


dq^ 


+ 


lj2=0 


cwsw 


dA^z 


dq'^ 


, (4.3) 


(j 2=0 


(4.4) 


Here, a = e^/ ( 47 r) = ( 47 r) is the fine-structure constant, Sw = sinS^ and = cos^w 

are the sine and cosine, respectively, of the Weinberg angle and the functions Ayv {q'^) 
are the coefficients of in the vacuum-polarization tensors Hy^, (q) = g^'^Ayy ((?^) + 
q^q'^Byyi (g^), where VV could be either 77 , 7 Z, ZZ, or WW. In our model, the oblique 
parameters are given by [ 22 ] 


AT 


’A [-f {«i. 0 - ("S'. <)] + 4 [F (Ml Ml) - F (Ml, Ml)] 

- [F (Ml, Mi) - F (Ml, Ml)] } , (4.5) 


m 2 


= > ^ 1 ) + {Ml ,Ml,Ml+clln^ + si In 




h2 




Mi 


+4 G (Ml , Ml) + si a (Ml ,Ml)-G (Ml Ml )} , 


(4.6) 


where the functions F, G and G are given in the appendix and M/j = 125.09 GeV denotes 
the reference value. 


4.3 Higgs invisible decay 

The model can also face constraints from the invisible Higgs decay, B{h —>■ inv) < 17% [23]. 
In our case we have inv = {/i 2 ^ 2 }) {-^dm-^dm}; when kinematically available. The corre¬ 
sponding decay widths are given by 

r (5, ^ ^ (^1 - e (M,. - 2Ms .). 

r {hi -4 N^mN^m) = 0 {Mh, - 2Mdm) . (4.7) 

The effective cubic coupling A 122 is defined below in Eq. (5.8). Due to the SI symmetry, we 
find that A 122 vanishes at tree-level, with the small (loop-level) coupling sufficient to ensure 
the decay to /12 pairs is highly suppressed.^ 

®Note that /12 decays to SM states, much like a light SM Higgs boson but with suppression from the 
mixing angle, s\. However, currently there are no dedicated ATLAS or CMS searches for light scalars in 
the channels 2 b, 2 r or 2 y, so we classify the decay hi h2h2 as invisible. In practice the suppression of 
r(/ii —>■ 112/12) due to SI symmetry renders this point moot. 



















4.4 The Higgs decay channel /i —)• 77 


The existence of extra charged scalars modifies the two Higgs branching ratios i3(/i —> 
77 , 7 Z), and this deviation can be parameterized by the ratios: 


_ B{h 77 ) 
^77 ^ 

" 7 ^ ^ 7 Z) 


^ 2Ch {tw) + A^cQf^i/2 

Sw'y {rs2,^S2) 

Ch (tw, Aw) + 


(4.8) 


(4.9) 


where tx = M|^/4M^ and Xx = M|/4M|^, with Mx is the mass of the charged particle 
X running in the loop, Nc = 3 is the color number, Qt is the electric charge of the top 
quark in unit of |e|, and the loop amplitudes Ai for spin 0 , spin 1/2 and spin 1 particle 
contribution [24], which are given in the appendix. Here 'dj, are the SM-like Higgs couplings 
to the pairs of charged scalars S^ 2 ) which are given by 


C/^AjjaT T SfiXfpQ^X. 


(4.10) 


The effect of the charged scalars on (4.8) and (4.9) depends on the masses for 5^, the 
sign and the strength of their couplings to the SM Higgs doublet and the neutral singlet 
and on the mixing angle 9h- One can use the reported results from LHC to constraints 
these parameters. 


5 Dark Matter 
5.1 Relic Density 

The lightest Z 2 -odd held is a stable DM candidate. As mentioned already, the lightest 
exotic fermion A^dm = is the only viable DM candidate in the model. The relic density 
is given by [25] 


^Duh 


2 


1.04 X lO^GeV^ 1 

Mpi ^/g*{Tf) < avr{xf) >’ 


(5.1) 


where Mpi = 1.22 X 10^® GeV is the Planck scale, g^{T) is the total effective number of 
relativistic particle at temperature T, and 


(cr(A/3M NY,-M)Vr) — 'y ^ (<T(iVoM A/jm —>■ X)Vr) 


X 




T ) 


[ ds CTNum (s y/sKi T j , (5.2) 


is the thermally averaged DM annihilation cross-section, Vr is the relative velocity, s is the 
Mandelstam variable, are the modihed Bessel functions and crjVoM XuM^allis) is the 
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annihilation cross into all kinematically accessible final state particles at the CM energy 

Vi. 


The parameter Xf 
puted from 


^DM/Tf represents the freeze-out temperature, and can be com- 


X/ =ln 


0 . 03 Mp/Mdm < (TVr{xf) > 


(5.3) 


As will be discussed in the next section, we require that to be in agreement with 

the observed value of the dark matter relic density [26]. 

The thermally averaged annihilation cross-section can be approximated in the non- 
relativistic limit as < aVr >= a -|- bv^, where Vr is the relative DM velocity and a and b 
are the s-wave and p-wave factors, which receives contributions from different annihilation 
channels. In this limit, the velocity squared is approximated by ~ 6/x/. Here, we 
evaluate the thermally averaged cross section exactly following (5.2). 


5.2 Annihilation cross section 

In our model, there are many contributions, where the channels can be classified into 
three types according to their Feynman diagrams types: (1) annihilation into charged 
leptons A^dmAom —t (Fig. 2-a and -b), which are t-channel diagrams mediated by 

charged scalars® (2) annihilation into SM fermions and gauge bosons pairs A^dmA^dm 
//, W~W'^, ZZ (Fig. 2-c), which occur through s-channel /ii^ 2 -mediated diagrams, and 
(3) the annihilations into scalars, A^mAom —t hi^ 2 hi ^2 (Fig. 2-d, -e and -f), which occur 
through both s- and t-channel diagrams. 



id) (e) (/) 


Figure 2. Different diagrams for DM annihilation. 

Charged leptons annihilation channel 

The DM Ni couples to SM leptons through the Yukawa couplings gia, and can annihi¬ 
late into charged lepton pairs as shown in Fig. 2-a and -b. The cross section for annihilation 

®Actually, for the same flavor case there two s-channel diagrams mediated by hi,2, however we neglect 
them due to the suppressed Higgs charged leptons couplings. 


10 - 


















into charged leptons ^ is given by [27] 


(-^DM-^DM 


8 

"^3 


1 \9lag*i0\^ 

87rs(M|+-M2^ + f)2 

“ -^dm)^ + ^ i ^ s + ~ -^dm) + % S /'S 
(M|+-M2, + f)2 4U ^ 


rrio 


+ ^L 




(5.4) 


SM fermions and gauge boson channels 

The processes A^dm^^dm bb, ti, W~^W~ and ZZ can occur as shown in Fig. 2-c. The 
corresponding amplitude can be written as 


M 


ichShyiu{k2)u{ki) 




Mh^SM {ruh Vs) , 


(5.5) 


where A4h^SM {fnh —^ \/s) is the amplitude of the Higgs decay h —)■ XsmXsm, with the 
Higgs mass replaced as —>■ V^. This leads to the cross section 


a{NuMNoM —t XsmXsmVt 


^V~sslclyl 


1 

»- K, 



2 

^h^XsMXsM Vs) , 

(5.6) 


where ^h^XsmXsm Vs) is the total decay width, with mh —>■ Vs- 

Higgs channel 

The DM can self-annihilate to /ii, 2 ^ 11 , 2 ) as seen in Fig. 2-d, -e and -f. The amplitude 
squared is given by 


\M\^ = ^fuMS 


^h^lik ^ ^h^2ik 


-\ 2 


»- K 


+4cj Cfc yjjjy[ i4f] 

{t-MU 

2 ‘^cf4yt 


^ - K 


^h^lik ^ ^h^2ik 


s - Ml 


0 - K 


',-Ml + Ml s + Ml-Ml 

- + a- ^ 




u-ML 


+ 


2 {4M2^m2^ + (m2^ + Ml - t) (m2^ + Ml -u)- sMl) 


+ (^dm + Ml - u) [Ml^ + Ml -t)- sMl ] 


t - MlV 

+ (m2^ + Ml - u) {Ml^ + Ml -u)-{s- 4M2^) {s - Ml - Ml)] , (5.7) 




with s, t and u being the Mandelstam variables, the Yukawa couplings defined as = yi, 
Cl = Ch and C 2 = Sh- Here, we integrate numerically on the phase space in order to get the 

^Indeed for same flavor charged leptons (a = / 3 ), there are hi,2 mediated s-channel processes that are 
proportional to their Yukawa couplings; we ignore these due to the Yukawa suppression. 
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cross section for a given s value. At tree-level the effective cubic scalar couplings (Aij^ and 
X 2 ik) are given by 

Am = 6 Ah civ - 3X^nclshV + SX^^Chslx - GX^slx, 

Aii2 = X^uclx + 2clsh{3Xu - X^a)v + ‘ichsl{2>X^ - X^^)x X^^slv, 

A222 = A122 = 0 , ( 5 . 8 ) 

though for completeness we use the full one-loop results that can be derived from the loop- 
corrected potential following [28]. The absence of cubic interactions and /z-l, at leading 
order, is a general feature of SI models. 


5.3 Direct Detection 


Concerning direct-detection experiments, the effective low-energy Lagrangian responsible 
for interactions between the DM and quarks is given by 

= ciq qq N^j^NY)yi, (5-9) 

with 


ShChMqXloM 

2vx 


K 


Ml 


(5,10) 


Consequently, the nucleon-DM effective interaction can be written as 




(5.11) 


with 

_ SfiCh {Mm — l-Alg) AfoM 


In this relation, Mj^ is the nucleon mass and Mg the baryon mass in the chiral limit [29]. 
Thus, the approximate expression of the spin-independent nucleon-DM elastic cross section 
at low momentum transfer reads 

(Mdm + Msf 

As will be discussed below, the most stringent constraint on a^et comes from the present 
as well the recent upped bound reported by LUX experiment [30, 31]. 


Ml 


Ml 


(5,13) 


Ml 


Ml 


(5,12) 


6 Numerical Analysis and Results 

In our numerical scan we enforce the minimization conditions, Eqs. (2.7) and (2.8), vacuum 
stability, the Higgs mass M/jj = 125.09 =F 0.21 GeV, as well as the constraints from LEP 
(OPAL) on a light Higgs [32]. The constraint from the Higgs invisible decay B{h inv) < 
17% [23] is also enforced. All dimensionless couplings are restricted to perturbative values 
and we consider the range 200 GeV < {(p) < 5 TeV for the beyond-SM VEV. We find a range 
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of viable values for consistent with the OPAL bounds, as shown in Figure 3. For the 
parameter space in our scan we tend to hnd in the range 0(1) GeV < < 90 GeV. 

Lighter values of appear to require a degree of engineered cancelation among the 
radiative mass-corrections from fermions and bosons; see Eq. (2.17). We noticed that 
regions with {(f)) > 700 GeV tend to be preferred in our scans. 



0 20 40 60 80 100 


Mh, (GeV) 


Figure 3. Scalar mixing versus the light scalar mass. The palette gives the branching ratio for 
invisible Higgs decays, with an overwhelming majority of the points shown satisfying the constraint 
B[hi —>■ inv) < 17%. 

We also scan for viable neutrino masses and mixing, subject to the LFV and muon 
anomalous magnetic moment constraints, while also demanding a viable DM relic density. 
In Figure 4 we plot viable benchmark points for the Yukawa couplings gia and /q,^, along 
with the corresponding LFV branching ratios and (5a^ contributions. It is clear that the 
couplings /a /3 are generally smaller than the couplings gta, and that the bound on r —)• /xy 
is readily satished, while the constraint from /x —)• eq is more severe. We observe in Figure 4 
that the model requires the largest coupling in the set gia to take 0(1) values. This feature 
is a generic expectation for three-loop models of neutrino mass, as one cannot make the 
new physics arbitrarily heavy, while reducing the Yukawa couplings, and retain viable SM 
neutrino masses. Thus, the testability of such models, which predict new physics at the 
TeV scale, is generally coupled with a need for 0(1) couplings. Consequently one expects 
such couplings to encounter a Landau pole in the UV, requiring a new description. We 
note that, when considering only one or two generations of singlet fermions, no solutions 
that simultaneously accommodate the neutrino mass and mixing data, low-energy flavor 
constraints, and the DM relic density, were found. Therefore at least three generations 
of exotic fermions are required. Also, we verified that the constraints from neutrino-less 
double-beta decay searches are easily satisfied for all benchmark points. 

Recall that, with regards to the DM relic density, there are many classes of annihilation 
channels, namely ^ X (X = bb, tt, WW + ZZ, /ii^ 2 ^i, 2 )- According to the 

DM mass, each channel could be significant or suppressed. In order to probe the role of 
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Figure 4. Left: Viable benchmark points for the Yukawa couplings gia and fap, in absolute values, 
where the dashed line represents the fully degenerate case, i.e, min |/| = max |/|. Right: The LFV 
branching ratios, scaled by the experimental bounds, versus the muon anomalous magnetic moment. 
The vertical line represents the muon anomalous magnetic moment experimental constraint. 


each channel, we plot the relative contribution of each channel to the total cross section, i.e. 
the ratio at the freeze-out versus the DM mass, in Figure 5-left. We see that the 

channel A^dm-^dm always fully dominant except for a few benchmark points. For 

DM masses smaller than 80 GeV the contribution oi X = hb can be significant, while in the 
range between 80 GeV< M^m < 100 GeV, both gauge bosons X = WW -|- ZZ and X = tt 
contributions can be important. In the range 200 GeV< M^m < 400 GeV, their contribution 
can reach 20%. For large DM masses M^m > 200 GeV, the X = hh contribution can reach 
at most 8%. The fact that the X = ti contribution could be important around 100 GeV, 
i.e., for Mom < M^, can be understood due to thermal fluctuations. Figure 5-right shows 
the corresponding charged scalar masses. For lighter DM masses of Mqm < 300 GeV, the 
charged scalar masses Ms ^2 should not exceed 450 GeV, while for larger values of 
the scalar masses Ms^ 2 ^e at the TeV scale. Such light charged scalars can be within 
reach of collider experiments [33]. 

Next we discuss the constraints from direct-detection experiments. We plot the direct- 
detection cross section versus the DM mass for our benchmark points in Figure 6. One 
observes immediately that the direct-detection limits impose serious constraints on the 
model, with a large number of the benchmarks excluded by LUX [30] as well the improved 
LUX bounds [31]. We find that only few benchmarks with Mom ^ 10 GeV or Mom ^ 
400 GeV survive the LUX bounds. As is clear from the figure, the surviving benchmarks 
will be subject to future tests in forthcoming direct-detection experiments. The palette in 
Figure 6 shows the corresponding values for M/ij, in units of GeV. In the region of parameter 
space for which Vom gives viable dark matter, we find that the Mh^ must be greater than 
20 GeV. 

We emphasize that we only found a few benchmarks for which the DM relic density 
was primarily determined by annihilations into scalars. On the surface, this claim may 
appear contrary to the results of Refs. [34, 35], which consider Majorana DM coupled to a 
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Figure 5. Left: The relative contributions of each channel to the annihilation cross section at 
the freeze-out temperature versus the DM mass. Right: The corresponding charged scalar masses 
versus the DM mass. 



1 10 100 1000 


Mdm (GeV) 


Figure 6. The direct detection cross section versus the DM mass compared to the recent results 
from LUX. The palette shows the mass for the neutral beyond-SM scalar, Mh^, in units of GeV. 

singlet scalar that communicates with the SM via the Higgs portal (called the Indirect Higgs 
Portal [34]). Naively one may expect our model to admit parameter space where the DM 
relic-density is determined primarily by the annihilations A^dm A^dm hh, in analogy with 
the results of Refs. [34, 35]. However, due to the SI symmetry, our model contains no bare 
mass terms, which reduces the number of free parameters in the Lagrangian. Consequently 
the DM mass M^m is related to both the coupling between Ndm and cj), and the mixing 
angle 9^. This reduction in parameters means we cannot evade the LUX constraints whilst 
generating a viable relic density by annihilations into scalars, explaining the difference 
between our results and Refs. [34, 35]. It also explains some features of the benchmark 
distributions in Figure 6. The benchmarks with larger contributions from the channel 
A^dmA^dm hh have a stronger coupling between A^dm and 4>- This increases the direct- 
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detection cross section due to hi ^2 exchange, creating conflict with the bounds from LUX, 
so the corresponding benchmarks are strongly ruled out. Indeed, with the smaller number 
of parameters in the SI model, it is a non-trivial result that viable regions of parameter 
space were found in Figure 6 . 




AS R(h->YZ) 

Figure 7. Left: The oblique parameters AS” versus AT for the benchmarks used previously. The 
palette shows the mixing sin^ 9h and all the points are inside the ellipsoid of 68 % CL. Right: Ratio 
of the widths for /i —>■ 77 and h 7 Z relative to the SM values. The constraints from ATLAS and 
CMS are shown, along with projected sensitivities after Run II at the LHC. 

Finally, we mention that the exotics in the model allow for new contributions to the 
Higgs decays /i ^ 77 and h —)• 'yZ. We plot the ratio of the corresponding widths relative 
to the SM values in Figure 7-right. We observe that a signihcant portion of the benchmarks 
are consistent with existing constraints from ATLAS and CMS. Importantly, the model can 
be probed through more precise measurements by ATLAS and CMS after Run 11. We note 
that all benchmark points are consistent with the oblique parameter constraints, as shown 
in Figure 7-left. 

7 Conclusion 

We presented a scale-invariant extension of the SM in which both the weak scale and 
neutrino mass were generated radiatively. The model contains a DM candidate, in the form 
of a sterile neutrino A^m- A new light neutral scalar is also predicted, namely the pseudo 
Goldstone-boson associated with the broken scale-invariance, / 12 , along with two charged 
scalars Si^ 2 - The masses for the latter are generically expected to be near the TeV scale, due 
to the related birth of the exotic scale and the weak scale via dimensional transmutation. 
The constraints on the model are rather strong, particularly the direct-detections constraints 
from LUX. However, we demonstrated the existence of viable parameter space with M^m ^ 
10 GeV or Mqm ^ 400 GeV. The model can be tested in a number of ways, including 
future direct-detection experiments, collider searches for the charged scalars, improved LFV 
searches, and precision measurements of the Higgs decay width to neutral gauge bosons. 
We note that the model does not possess an obvious mechanism for baryogenesis - it would 
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be interesting to study this matter further. In a partner paper we shall study the scale- 
invariant implementation of the Ma model in Ref. [4]. 
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A Multi-Scalar Scale-invariant Theories 

In a general multi-scalar theory one cannot minimize the full one-loop corrected potential 
analytically. However, with recourse to the underlying SI symmetry, there exists a simple 
analytic approximation that captures the leading features [15]. A general tree-level SI 
potential for a set of scalars {(Pa} can be written as 

kb({</’A}) = gABCD (pA(pB (pc (pD, (A-1) 

where the dimensionless couplings gABCD are symmetric. In general, these couplings are 
running parameters that depend on the energy scale, gABCD = gABCD{g)^ and one can 
freely select a value of jj, that simplifies the analysis. A convenient choice is the value 
// = A, at which the tree-level potential vanishes along the direction of an assumed non¬ 
trivial minimum in field space, namely 

gABCD (A) (pA (pB (pc pD = 0. (A.2) 

Here, the minimum is defined by (Pa) = RpA, with pA a unit vector in field space and R 
a (yet to be determined) radius. Combining Eq. (A.2) with the minimization conditions, 
dVo/dpA = 0, determines the angular VEVs pA in terms of the couplings gABCD- Subse¬ 
quently expanding around the ground state in the tree-level potential reveals a spectrum 
containing a massless scalar, corresponding to the flat direction. 

Eq. (A.2) implies that the tree-level potential vanishes, at the scale /i = A, to an 
accuracy on the order of the loop corrections: 

V^{{RPA}]g = A) < 0(Ri_ioop), (A.3) 

where we display the renormalization scale dependence and write the full loop-corrected 
potential as R = Vq + ki_ioop + ■ ■ ■■ Thus, one-loop corrections can be comparable to Vq 
along the direction pA, so the interplay of the two terms allows a non-trivial minimum that 
lifts the flat direction to fix the radial VEV (R). Adding the one-loop corrections along the 
direction pA gives 

V{{RPa}; g = A) = VoURPa}-, g = A) + Vi_ioop{{RpA }\= A) + ..., (A.4) 
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which can be written as [15] 




d 2 


(A.5) 


with 


A = 


64:Tt'^{R)^ 


Tr 


Ms log 


{Ry\ 


-Tr 


, >4 1 -M^ 


+ Tv 


A ^4 1 -^V 

My log 




, (A. 6 ) 


and 


Here A4s,F,v th® mass matrices for scalars, fermions and vectors, respectively, and 
the trace runs over both particle species and internal degrees of freedom. Minimizing the 
one-loop corrected potential lifts the flat-direction to give 

(E) = g-(A/2B+l/4}j^^ 

The dilaton acquires a loop-level mass, given by = 8I3{R)‘^. Thus, radiative correc¬ 

tions successfully induce a non-trivial VEV for one or more of the scalars (j)A, by introducing 
a dimensionful parameter, (R) oc A, in exchange for one of the dimensionless couplings in 
Eq. (A.2). This manifests dimensional transmutation. 

In the present model, demanding that = 8B{R)‘^ > 0, requires that B be 

dominated by the term TrAdg, meaning that one (or both) of the scalars 3^2 must be the 
heaviest state in the spectrum. In practise, this implies that A is also dominated by the 
contribution of 5]’’2 to the Adg term in Eq. (A. 6 ). Thus, loop corrections from the scalars 
hi ^2 along the flat direction are sub-dominant to the corrections from 3^2-^ Therefore, 
simply dropping the corrections from hi ^2 will not introduce a significant error in the analysis 
(the error is expected to be 0{M^^/Ms^ ^)). As discussed in the text, this simplification 
has the advantage of allowing one to obtain analytic expressions for the ground state by 
minimizing the one-loop corrected potential directly. As a point of comparison, for the 
present model, the minimization in Eq. (A.2) gives Ah (A) A 0 (A) -|- A,^h(A) = 0, and we 
see from Eq. (2.7) that our approach incorporates loop corrections to this expression, up to 
0{M^_^/Ms-^ ^) effects. Taking the heaviest scalar as Ms > 300 GeV (which we can always 
do - see Figure 5), the error in the loop terms is typically < 3%. Once we have found the 
ground state, we reintroduce loop corrections from hi ^2 to determine the mass eigenvalues, 
reducing the error in the expressions for the scalar masses and mixings. 


®For parameter space of interest in this work, corrections from are also smaller than those from the 
top quark and, in large regions of parameter space, one or more of the fermions N. 
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B Oblique Parameter Functions 


The functions employed in the calculation of the oblique parameters in Section 4 are defined 
as follows: 


F{I,J) = 


1+1 _ LL In -t <^= f ^ 7 

2 I-J j ^ J ^ 
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^ I = J, 


(B.l) 
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with t = I + J — Q and r = — 2Q (/ + J) + (/ — J)^ , and 
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y/r In 
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r > 0 , 
r = 0 , 
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_ 2 arctan 

C Loop induced Higgs decay functions 

The functions used to evaluate the Higgs decay rate of h —)• 77 are given by 

^ 0 ^ {x) = -x~‘^ [x- f (x)], 

^ 1/2 (®) = [x + {x-l)f (x)] , 

Aip (x) = —x“^ [2x^ + 3x + 3 (2x — 1) / (x)] , 


with 


fix) = 


f 

arcsin^ (\/^) 

IH 

log in 


X < 1 
X > 1, 


and those used in the decay rate of /i —)• 7 Z are given by 

^ 0 ^ ix,y) = h ix,y ), 


with 


(B.3) 


(B.4) 


(C.l) 


(C.2) 


1/2 (x, y) = h ( 3 :, y) - h (x, y ), 

Al^ (x, y) = [(1 + 2x) tan^ 0^ - (5 + 2x)] h (x, y) + 4 (3 - tan^ 6^) h (x, y) , (C.3) 

T — 1 I lA-fiy) I yl9A-g{y)] j _ fA-f(y) (c 

H y) - + 2 {x-yr + {x-yf ’ (3:, yj - 2 (,:_y) , l<-4j 
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and 
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\/X ^ — 1 arcsin {\fx) x < 1 



(C.5) 
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